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ABSTRACT 


A. new approach to the eigensystem assignment problem is presented. 

The approach utilizes a null-space formulation of the eigenvalue/eigenvector 
assignment problem to simultaneously realize arbitrary eigenvalue specifi- 
cations, approximate desired modal behavior, and achieve low eigensystem 
sensitivity with respect to plcint parameter variations. The methods are 
applied to the design of regulator and integral plus proportional servo 
control systems. 
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A SUMMARY OF SPECTRAL SYNTHESIS PROCEDURES 


FOR MULTIVARIABLE SYSTEMS 

By 

S.R. Liberty^, R.R. Mielke^, and R.A. Maynard^ 

DISCUSSION 

State space-oriented procedures provide significant design advantages 
over classical control design methods for flight control systems. Several 
procedures involving the synthesis of multivariable systems based on eigen- 
value-placement criteria have previously been developed (refs. 1 to 3) . 
Eigenvalue specifications alone do not characterize the actual variable 
response of a system, however. The eigenvector corresponding to a given 
eigenvalue determines the influence of that eigenvalue on the state variable 
response (ref. 4), Thus, satisfying modal matrix specifications is as 
important as achieving eigenvalue specifications in designing a feedback 
system. 

The multivariable design problem can therefore be approached as an 
eigenvalue/eigenvector assignment problem (refs. 5, 6). This approach has 
led to a null-space formulation of the problem (ref. 7) . All possible 
modal matrices related to a given eigenvalue specification have been shown 
to lie in the null space of a corresponding matrix. A desired modal matrix 
is projected onto this null space, yielding the closest approximation to 
the desired modal matrix that is actually realizable. 

System sensitivity to plant parameter variations is then reduced 
local to the modal matrix obtained from this projection (ref. 8) . This is 
done by perfo.rming a gradient search within the null space for a modal 

^Professor and Chairperson, Department of Electrical Engineering, Old 
Dominion University, Norfolk, Virginia 23508 

^Associate Professor, Department of Electrical Engineering, Old Dominion 
University, Norfolk, Virginia 23508 

^Graduate Research Assistant, Department of Electrical Engineering, Old 
Dominion University, Norfolk, Virginia 23508 


matrix that minimizes a sensitivity cost function and that lies near the 
modal matrix first obtained. 

The null-space formulation of the eigenvalue/eigenvector assignment 
problem is applicable to the design of regulator and servo control systems 
(refs. 9 and 10). Dynamical compensators or observers (ref. 11) are used 
to realize the integral plus proportional feedback law determined under the 
assumption of complete state accessibility. The final system is such that 
the plant/integrator eigensystem and observer eigensystem are decoupled and 
can be independently specified. 

Computer programs based on the above procedures have been written and 
are currently being evaluated. Initial results indicate that control 
systems having desirable modal behavior and low eigensystem sensitivity to 
changes in plant parameters can be designed using these methods. A diffi- 
culty with this procedure, and all other eigensystem assignment procedures, 
is the uncertainty of how to translate system performance specifications 
into specification of a desired modal matrix. Continued research in this 
area would appear to be fruitful. A complete and detailed development of 
all procedures mentioned here is currently being prepared for submission 
to a refereed journal (ref. 12). 
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ABSTRACT 

A method for selecting a multivariable state feedback controller that 
simultaneously achieves an a priori specification on closed loop eigenvalues 
and good mode mixing is presented. The problem Is solved by projecting a 
desired modal matrix onto a constraint set containing the null space of the 
closed-loop state matrix, while assuring that the projection 1« In the null 
space. The feedback matrix follows Inroediately in the formulation. An 
example Involving a helicopter hover controller Is presented, 

I. INTRODUCTION 

This paper Is concerned with the design problem of selecting the 
stationary feedback matrix K in the control law u » Kx for a controllable 
multivariable system x » Ax + Bu where xcR" and ueR"*. The matrix K Is to 
be selected so that the closed-loop system matrix A + BK has arbitrarily 
assigned eigenvalues. Once the eigenvalues are assigned, the remaining 
freeoom available In selecting eigenvectors is utilized to achieve desir- 
able modal behavior. Srinathkumar and Rhoten In [1], [2], and [3] have 
presented a "spectral synthesis algorithm" for constructing the feedback 
gain K which exploits the nonuniqueness of the modal matrix of A + BK In an 
attempt to achieve a priori design specifications. However, the a priori 
specifications are not explicitly incorporated in the procedure and th« 
Inherent design freedom available is not systematically exploited. 

The work reported here relies heavily on Insight gained in [1] but is 
based upon a different approach to the synthesis problem that allows a 
priori modal specifications, and results in a simple design procedure that 
can be parameterized and iterated, 

II. NULL SPACE FORMULATION OF THE PROBLEM 

The eigenvalue/eigenvector assignment problem can be expressed In the 
form AU ■ UA, where A * A + 3K, U Is the modal matrix and A is a diagonal 
matrix whose nonzero entries are the desired eigenvalues of A. 


This work was supported by the NASA-Langley Research Center under 
grant number NSG-1519. 
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This ctn bt expressed In the form 

CdnOA) + {-A^« !„)] vec U • 0, .0) 

where® Is the Kronecker product [4] and vec U ■ col (uj). Thus, If IJ Is a 
solution to • UA, then vec U lies In the null space of S • C(I„0A) 

*■ Moreover, S • vec U ■ 0 If and only If [A - U| " 0 

for all 1 ■ 1 , •••» n. 

In the remainder of thli^i paper It is assumed that the system is in 
rank-reduced form with B « Cljj, | 0]^, If this is not the case, a coordinate 
transformation is applied to the system so that it is in this form. A 
solution is obtained according to the procedure to be described, and this 
solution is then transformed back to the original system coordinates. If A 
Is partitioned so that 


*11 

*12 

*21 

* 22 . 


with A^^ an m X m matrix, then [A - X^Ij^] u^ ■ 0 if and only if the following 
two equations hold; 

(CA„ - X,IJA^2^ + K) u, - 0 ' (2) 

and 

[Aj^lAjg - u^ - 0 (3) 

Equation (3) is a constraint on the vectors that lie in the null space of 
[A - X^Ijj]; that is, vectors in the null space of [A - X^I^] must also lie 
in the null space of CA 21 IA 22 - Equation (2) is regarded as a 

constraint on K to be satisfied after U has been computed. 

The desired modal mix is expressed as a desired modal matrix, P. By 
the projection theorem [5], the closest approximation, in a least squares 
sense, to P that can be achieved is obtained by projecting vec P onto the 
null space of S. This is equivalent to projecting p^ onto the null space 
of [A - X^I^] for each 1 » 1 , ••• , n to get each of the columns of the 
modal matrix, U, that most closely approximates P. This projection pro- 
cedure can result in a U that is not a valid modal matrix. This occur- 
rence will be indicated by the noninvertibility of U. When this happens, 
the entries of P are perturbed and the new P is projected. Once an in- 

A A 0 1 A 

vertible U is obtained, A is computed as A = UAU with A partitioned to 
conform with A. The feedback matrix K is calculated from 



Coffiplex eigenvalues are treated In essentially the same way as real 
eigenvalues, after some preliminary manipulation. Using ± J8^ 

end • X, * jy^, Au, • u^A^ becomes A(Xj t jy^) • (x, t jy^)(a^ * jp^). 
The real and Imaginary parts of this equation can be separated and the 
resulting equations can be written as 


Partitioning 


r » « T 
“1 *1 

AC*,ly^3 ■ [x,|yj 

L '®1 “ 1 . 



-(5) 

( 6 ) 


where s^ and t^ represent the first m-components of x^ and y^, respective- 
ly, we obtain two constraint equations on u^! 


Aji s^ i A 22 - c.^ v^ - w^ , (7.) 

*21 *1 *22 ^1 *8 ^ V ^ + w ^ . ( 8 ) 

These can be expressed In the form 



Thus, the constraint on u^ and Its conjugate Is that the 2-n component vec- 
tor (x^ yp' must lie In the null space of 



Letting p^ and q^ denote the real and Imaginary parts, respectively, of the 
desired eigenvector conjugate pair associated with s js^, (pj q|)^ Is 
projected onto the null space of to obtain the columns for the real and 
Imaginary parts of the corresponding eigenvector conjugate pair In the sys- 
•tem modal matrix. 
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m, HELICOPTER HOVER COHTROLLER APPLlCAflOr* 


The design procedure 1$ illustrated by an application te a problem 
considered by Srinathkumar [I3i The Sikorsky SH-30 helicopter In a hover 
mode is modeled as a ninth order linear system having state variables of 
longitudinal (u), lateral (v), and vertical (w) velocities in feet/seconds*, 
pitch (q), roll (p) and yaw (r) rates in degrees/seconds; and pitch (s), 
roll (♦) and yaw (o) angles in degrees. The Inputs are main "'tor collec- 
tive pitch (Uq), tail rotor collective pitch (u^), longitudinal cyclic 
pitch (Up) and lateral cyclic pitch (u,^), all in degrees. 

The open-loop dynamics are given by 

^1 , ^11 *12 ^1 ^ ‘^11 ®12 “l 

[ij] Aj, Ag2 Xj Bg, B22 Ug 

where ■ [o ** R *nd Ui • [Up u^]^ are the longitudinal variables and 
controls, and X 2 * [v p a r *3^ and Ug « [uj^, Uj]”^ are the lateral variables 
and controls. The normalized j/stem matrices, scaled by a rotor tip speed 
of 680 ft/sec, are gi^en in Table 1/ along with the open-loop eigenvalues. 
The desired slosed-lgop eigenvalues are: 

Response Variable Eigenvalue 


U 

w 

q.e 

V 


p.'/ 

r,* 


-4.5 
-.324 
-1.5 i J1 
-.3 

-1.5 t jl 
-1.5 ± jl 


With mode decoupling as an objective, a block diagonal P matrix is 
specified. The resulting modal matrix U and control law matrix K are given 
1n Table II. As a second example, another block diagonal P matrix and the 
resulting U and K matrices are given in Table III. 


If the modal behavior resulting from the above procedure Is unsatis- 
factory, then utilization of a weighted norm in the projection procedure 
allows the designer to Iterate to a more satisfactory solution for the pre- 


selected P matrix. For example, if an element of the modal matrix resulting 
from a projection is too large, then the norm over which optimization is 
carried out (via the projection theorem) should be weighted more heavily in 
the corresponding component. This weighted norm idea, which is not exem- 
plified here, allows one to carry out an iterated design procedure. 
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TABLE II. FIRST ITERATION 


- 0.5 - 0.5 
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ABSTRACT 

A new method of selecting a multivariable 
state feedback controller Is presented. The 
resulting controller simultaneously realizes 
arbitrary closed-loop'^ngenvalue*, appro.'timates 
specified modal behavior and achieves low 
elgensystem sensitivity with respect to plant 
parameter variations. The method cnaracterizes 
a vector space slightly larger than the null 
space of the closed-loop system matrix and pro- 
jects a desired modal matrix onto this space. 
Sensitivity of eigenvalues and eigenvectors is 
then minimized local to the desired modal 
matrix using a gradient search technique. A 
tutorial example to illustrate the design pro- 
cedure is given. 


I. INTRODUCTION 

Consider the controllable multivariable sys- 
tem x ■ Ax + Bu where x e R^ and u e R*”. In this 
paper we are concerned with the eigenvalue/eigen- 
rector assignment problem of selecting a station- 
ary feedback matrix K in the control law u = Kx 
s.uch that tlie eigenvalues of the closed-loop sys- 
tem matrix A * A + BK are arbitrarily placed. 

^ce the eigenvalues are assigned, the remaining 
freedom available in selecting the closed-loop 
tystem eigenvectors is to be used to approximate 
tpecified modal behavior and to achieve inse.isi- 
tivity to plant parameter variations. 

A "spectral synthesis algorithm" for con- 
structing a feedback gain matrix which utilizes 
the nonuniqueness of the closed-loop system matrix 
an attempt to achieve specified mocal behavior 
“as presented in [1], [2], and [3]. It was also 
shown how sensitivities of eigenvalues and eigen- 
vectors could be calculated for a ouen design. 
>*»“ever, a priori modal and sensitivi'.y specifica- 
tions were not explicitly incorporated in the pro- 
tadure and the inherent design freeccm available 


was not fully exploited. In [4'J the total design 
freedom available In selecting the eigenvectcrs’of 
A was used to globally minimize elgensystem sensi- 
tivity. However, the procedure often results in 
unsatisfactory system dynamic behavior since real- 
ization of low system eigensensitivity and desired 
modal behavior are sometimes conflicting objec- 
tives. 

The work reported here relies heavily on In- 
sight gained in [1] but is based upon a different 
approach to the synthesis problem. The first step 
inthe design procedure is the calculation of the 
initial closed-loop system eigenvector matrix 
which most closely approximates the modal matrix 
design specification. This is accomplished by 
constructing an orthonormal basis of a vector 
space slightly larger than the null space of the 
closed-loop system matrix and then projecting the 
desired modal matrix onto this space. The second 
step of the design procedure is to reduce eigen- 
system sensitivity to plant parameter variations. 

A cost function consisting of a weighted sum of 
eigenvalue and eigenvector sensitivities is mini- 
mized using a gradient search procedure. The 
search is performed local to the initial eigen- 
vector selection to retain desired modal char- 
acteristics. 


II. INITIAL EIGENVECTOR SELECTION 

It is assumed that the system is in rank- 

reduced form with B = [1^^^ | 0]^ and only the case 

of distinct eigenvalues is considered. The eigen- 
system equation for the closed-loop system is 

AU - UA (I) 

where U » [u^ Ug ‘ ' ‘ Uj^] is the modal matrix and 

A » dia {x^, Xg, ■ ' ', X|^} is a diagonal matrix 

of eigenvalues. It has been shown [5] that nec- 
essary and sufficient conditions for U to be a 
modal matrix satisfying (1) are that each eigen- 
vector u^ be in the null space of the matrix 


This work was supported by the NASA-Langley Research Center under grant number NSG-1519. 
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$1 • [Ajt 1 A 22 - X, wher« 


^1 

• 

Ai2 




CM 
f 

1 

1 

A 22 ^ 


and, that these eigenvectors be linearly indepen- 
dent in r”. Let be a matrix whose columns are 
an orthonormal basis for the null space of S^. 

Then the columns of all possible modal matrices 
satisfying (1) can be expressed as linear combina- 
tions of the columns of all R^'s. 

Let P denote the desired, although not nec- 
essarily realizable, modal matrix given as a de- 
sign specification. The initial eigenvector 
selection is made by projecting P, column by 
column, onto the vector space containing all pos- 
sible modal matrices. This results in n sets of 

constants, “21* “m1^^’ 

i » 1, 2, •••, n such that u^ ' R^ If the 
resulting eigenvectors U| are linearly indepen- 
dent, then the realizable modal matrix U which is 
closest to the desired modal matrix P in a Eu- 
clidean norm sense has been obtained. If the 
eigenvectors are not linearly independent or if 
the resulting U matrix is unsatisfactory, the 

parameters are perturbed slightly by v/eighting the 
projection norm. Once an invertible modal matrix 
is obtained, the closed-loop system matrix is com- 
puted as A = UaU*\ The initial feedback gain 
matrix K is then calculated from 

K “ I ^12^ * ^^11 ^ ^12^ , 

where A is partitioned to conform to A, 


III. EIGENSYSTEM SENSITIVITY REDUCTION 

At the completion of the first stage of the 
design process an initial feedback gain matrix K 
has been calculated such that the resulting 
closed-loop system realizes the specified eigen- 
values and approximates the desired modal matrix 
as closely as possible. In addition, the design 
freedom available in selecting a modal matrix U 
has been parameterized in terms of the constants 
a^.j. Suppose now that the elements of the plant 

matrices A and B are dependent on a scalar para- 
meter p. Then the closed-loop system is also 
dependent on p and the extent of this dependency, 
that is, the sensitivity of the system eigenvalues 
and eigenvectors with respect to p, is related to 
the choice of modal matrix U. Inthe second stage 
of the design procedure, the initial selection of 
U is modified- such that a cost function consisting 
of a weighted sum of eigenvalue and eigenvector 
sensitivities is reduced. 


Let ~ and ^ denote the sensitivities of 
plant parameters with respect to p, Then 


dp 


dA . dB 
dp * 2p 


K 


and the eigenvalue sensitivities are given 


( 3 ) 
by [ij 



( 4 ) 


and the eigenvector sensitivities are 



J-1 

m 


-^i ■ ^j _ 


( 5 ) 


where V » [v^ Vg ••• v^^] is a normalized recipro- 
cal basis for A. The sensitivity cost function Is 
defined by 



where -E., m^., i = 1 , 2, •••, n are positive weight- 
ing constants. Modifications to the initial modal 
matrix are calculated using a gradient search pro- 
cedure in which the sensitivity cost function J is 
reduced by choice of the constants a^j. Thus, at 

the k+1 iteration, a gradient 


TToferr 




(7) 

is calculated and used to generate new constant 
values, a^.j(k+l), from the current values, Ojj(k), 

according to 

aij(k+l) - a,.j(k) - (8) 


where d is the gradient search step size. 


Calculation of required derivatives of the 
cost function proceeds as follows. Defining two 
auxiliary sensitivity matrices S » [s^j] where 

S»U-Vf)U (9) 
and H ' [h.] hg ••• h^] « Cb^j] where 
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0 


h 


U 


_!iL 


1-J 


( 10 ) 


IV. A DESIGN EXAMPLE 

Consider the system x ■ Ax + Bu with x c 
ind u c R where 


the cost function Is expressed as 
n 


k-1 t J 


Straightforward differentiation with respect to a 
single coefficient then yields 


3J 


3a 


1 J 


■ 2 


where 


k-1 


i^k ®kk' 




( 12 ) 




!!k 

®* 1 j 


and 


3S 


VS \ «ij 


\ I 


3a 


U 


q-k 

q>*k 

.,/dBw3U^ 


(13) 




( 14 ) 


"2.0000 

0.0000 

1.0000 

1 

0.0000 

-(2.0000)/p 

(1.0000)/pj , 

_ 1.0000 

1,0000 

-2.0000' 

J 

"1.0000 

0.0000 



0.0000 

(1.0000)/p 



. 0.0000 

0.0000 _ 




and p Is a plant parameter with nominal value 
1.0000. It Is desired to calculate a feedback 
gain matrix K such that the closed-loop system has 
eigenvalues 

• -1.0000, X 2 - -1.2000, X 3 » -3.0000 
and modal matrix as close as possible to 


P 


3.7500 -0.6700 I.OOO 0 I 

3.2500 0.7500 -1.0000 : . 

7.0000 0.0000 O.lOOOj 


The first step in the design procedure Is 
the selection of an initial set of closed-loop 
eigenvectors. The eigenvectors are selected 
such that the initial modal matrix U is as close 
to P as possible without regard for elgensystem 
sensitivity. The vector space containing all 
possible modal matrices is characterized by the 
matrices 


0.7071 0.408f] P0,7071 0.34Sf| 
0.7071 0.4082 I . Rg = I 0.7071 0.3482! 
0.0000 O.8I65J j_0.0000 0.8704J 


The sensitivity reduction stage of the design 
process is highly designer interactive. At the 
completion of each iteration, a decision to con- 
tinue or to terminate the sensitivity reduction 
must be made. This decision often involves a 
compromise between achieving low eigensystem 
sensitivity and retaining desired modal behavior. 
Once the decision to terminate has been made, the 
final controller design is calculated using equa- 
tion ( 2 ). 

A computer program based on this design pro- 
cedure has been written and is currently being 
evaluated. Initial results indicate that con- 
trollers with significantly lower eigensystem 
sensitivity and satisfactory modal behavior are 
obtainable with this procedure. A tutorial ex- 
ample is presented in the next section to illus- 
trate the new design procedure. 


and 

p0.7071 -O'.AOsfj 
R3 0.7071 -0.4082 ' . 

[_0.0000 0.8165J 

Projecting P into this space yields the sefof 
a coefficients 


r -0.3536 

a( 0 ) - 

8.5732 

1.0041 

0.0279 

-1.4142*" 

0,0816 

the corresponding modal 

matrix 


“’3.7500 

-0.7003 

0.9667“’ 

U ■ 3.2500 

0.7197 

-1.0333 

j_7.0000 

0.0242 

0.0667 
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and tha Initial feedback gain matrix 


K • 


13.2526 

-13.1593 


12.5341 

-12.4526 


-13.3833 

12.2955 


For purposes of comparison with other modal 
matrices, U Is normalized so that 




0.4369 

0.3789 

0.8158 


-0.6972 

0.7165 

0.0242 


0.6324 

0.7295 

0.0471 


For this Initial modal matrix selection, the 
eigenvalue sensitivities are 

dX, ^ dx- ^ dXg ^ 

-0.11,|g^| “65.14, Ijpil - 492.69, 
and the eigenvector sensitivities are 


du, ^ du- ^ 

gplll - 18524.93, I Ig^l) - 39.75, 



With unity weighting coefficients for all eigen- 
value and eigenvector sensitivities, the total 
system sensitivity cost function is 

J » 19386.4, 

If the’ plant parameter p is perturbed lOS above 
nominal, then the closed-loop eigenvalues become 


V(l) 


0.0172 0.0040 -0.0578 

0.0094 -0.1372 -0.9887 


Using a step size d « 0.9000 results In a new set 
of a coefficients 


aO) 


•0.3690 1,0005 -1.3622 

8.5647 0.1514 0.9714 ’• 


a new modal matrix 


0.4383 

-0.6471 

0.3387 

0.3774 

0.7512 

-0.8127 

0.8157 

0.1302 

0.4741_ 

spending 

feedback 

gain matrix 

1.4904 

0.6961 

-1.5856 

-1.4195 

-0.6904 

0.5448 


For this selection of modal matrix, the eigen- 
value sensitivities are 


dX,^ dX, ^ dX, ^ 

la^l * 0.06, \-^\ » 0.01, Ig^l - 6.57, 

the eigenvector sensitivities are 



and the system sensitivity cost function Is 


X^ » -0.9915, Xg “ -1.4473 - j1.0981 , 
Xg - -1.4473 + J1.0981 


J ■ 19.2. 

If the plant parameter p is again perturbed lOS 
above nominal, the closed-loop eigenvalues are 


and the normalized modal matrix is 


X^ - -0.9767, Xg * -1.2104, Xg » -2.7683, 


0.6601 + jO.OOOO 
0.0918 + jO.OOOO 
0.7456 + jO.OOOO 


0,0000 jO.7262- 
0.0515 - jO.6834 
-0.0123 + j0.0531 


0.7262 + jO.OOOO 
-0.6834 + j0.0515 , 

0.0531 - jO.012,',1 


These results Indicate high eigensystem sensi- 
tivity to changes in the parameter p for the 
Initial eigenvector selection. 

The second stage of the design process con- 
sists of reducing the eigensystem sensitivity 
using a gradient search procedure. The gradient 
at the first iteration is 


and the normalized modal matrix is 




0.4479 

0.3802 

0.8092 


-0.7404 

0.6654 

-0.0949 


0.3717 

-0.7701 

0.5184 


The new results indicate significantly reduced 
eigensystem sensitivity has been achieved with 
relatively small changes in the Initial modal 
matrix selection. Continued iteration results 
in a final cost function value of approximately 
J * 8.65 with a correspondingly larger distance 
between the resulting modal matrix and the 
specified matrix P. 
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Abatract 

This briaf paper summarizes the approach the authors will take in designing a feed- 
back controller for the F-lOO turbofan engine. The technique to be utilized simul- 
eaneously realizes dominant closed- loop eigenvalues, approximates specified modal 
behavior, and achieves low eigensystem sensitivity with respect to certain plant 
parameter variations. 


SUMMARY 

Our approach to the design of a feed- 
back controller for the F-lOO turbo- 
fan engine [1] is essentially a pole- 
zero assignment technique. The actual 
controllers that we will design are in 
the multivariable proportional-plus- 
integral class. The matrix gains in 
those structures will be determined 
by manipulations on the eigensystem 

AU - UA 

where U is the modal matrix of A and 
A is the corresponding quasi-diagonal 
matrix of complex eigenvalues. The 
matrix A is a composite of feedback 
gains and original plant-model para- 
meters. Consequently, we call our 
technique a "spectral synthesis pro- 
cedure" or an "eigenvalue-eigenvector 
assignment technique". 

To date, this technique has only been 
applied in state feedback contexts 
where all closed- loop eigenvalues 
could be arbitrarily assigned. This, 
of course, is not the case in the 
P-lOO turbofan engine. There are two 
ways to handle this without modifying 
the existing state feedback technique. 
(Modification to the output feedback 
case is nontrivial and a subject of 
current research.) 

First, design could be carried out as 
if all states were available for mea- 
surement, and the resulting controller 
then coupled to an observer. Second, 
if the plant measurement, explicitly 
includes all dominant mode behavior of 
the system and nondominant modes are 
stable, then application of our tech- 
nique to the reduced order system where 
order equals the number of measurements 
is straightforward. 


The unique aspects of our eigenvalue- 
eigenvector assignment technique are 
the explicit incorporation of a pKioxl 
"mode mixing" specifications, and the 
systematic achievement of low eigen- 
system sensitivity within the mode mix- 
ing constraints. 

Low eigensystem sensitivity and desir- 
able mode mixing are often, but not 
always, competing objectives. The 
spectral synthesis technique, includ- 
i.ng these unique aspects, is documen- 
ted in [2] and [3], where it appears 
in the regulator context. 

Our oral presentation will illustrate 
eigensystem sensitivity for several 
designs and will contain full-scale 
simulations of these designs. Per- 
formance comparisons will be made and 
the design procedure will be demon- 
strated. 
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ABSTRACT 


II. THE EIGENSYSTEM ASSIGNMENT PROBLEM 


A new design procedure for a class of multi- 
variable integral plus proportional servo control sys- 
tems is presented, The procedure utilizes a null space 
formulation of the eigensystem as'iignftent problem 
assuming complete state feedback. Reduced order obser- 
vers are used to estimate unavailable states, The 
resulting system realizes arbitrary closed-loop eigen- 
values, approximates specified modal behavior and 
achieves low eigensystem sensitivity with respect to 
plant parameter variations, A tutorial example to 
illustrate the design procedure is included. 


I. INTRODUCTION 


A new design procedure for a class of multi- 
variable Integral plus proportional servo control sys- 
tems is presented, The design procedure consists of 
three steps. In the first step. Integral and propor- 
tional feedback gain matrices are calculated assuming 
complete state accessibility using an eigenvalue/eigen- 
vector technique. The technique utilizes a null space 
formulation of the eigensystem assignment problem which 
simultaneously realizes arbitrary closed-loop eigen- 
values, approximates specified modal behavior and 
achieves low eigensystem sensitivity with respect to 
plant parameter variations. The second step in the 
design procedure is the construction of a reduced order 
observer to estimate those states which are not acces- 
sible at the output of the plant. It is shown that the 
nuTl space formulation of the eigensystem assignment 
problem is also applicable in the design of reduced 
order observers, The final step in the servo problem 
design procedure is the calculation of system gain 
matrices to realize an ovc'.'ail eigensystem in which the 
plant/integrator dynamics and observer dynamics are de- 
coupled, 

The null space formulation of the eigensystem 
assignment problem is outlined in Section II. A more 
detailed presentation of the eigenvaluc/eigenvactor 
assignment technique is given in [1] and the method of 
reducing eigensystem sensitivity with respect to plant 
parameter variations is described in [<'.]. In Section 
III, the new design procedure for multivariable integral 
plus proportional servo control system!, is presented. 
This is followed by a tutorial example in Section IV to 
illustrate the design procedure. The interested reader 
is referred to [3] for a more detailed development of 
all procedures outlined here. 


Consider the multivariable system described by 

X » Ax + Bu whore x c R", u t r"', B has full column 
rank, and the pair (A,B) is controllable. The eigen- 
system assignment problem consists of solecting a 
real -valued mxn matrix K such that the closed-loop sys- 
tem matrix 


A - A + BK . (1) 

has arbitrarily specified eigenvalues and eigenvectors 
that are close in some sense to a desired eigenvector 
assignment. In addlticn, it is highly desirable that 
the eigensystem of A be insensitive to changes in ele- 
ments of A and B. In general it is not possible to 
achieve arbitrary eigenvector assignment so that it is 
necessary to seek an approximation to a desired assign- 
ment, The controllability assumptions on the pair 
(A,B) are sufficient to guarantee arbitrary eigenvalue 
assignmer.t. 

For clarity of presentation in this section it is 
assumed that a change of coordinates has been made such 
that B is in rank reduced form, B ■ [I^ 0], where 1^, 

denotes the mxm identity matrix. The eigensystem equa» 
tion for A is 


AU • UA 


( 2 ) 


where U » [u^ , Ug, •••, Uj,] is a matrix of eigenvectors 

of A and A ■ diagCXp Ag, •••, Aj,} is the diagonal 

matrix of eigenvalues of A. Only the case of distinct 
eigenvalues is treated here. If A and K are parti- 
tioned compatibly with that of B, 


A ■ 


.*21 I * 22 - 


K » [K^ I Kg] , 


then equation (2) is equivalent to 


and 


CA,i 


- A^ I + K,i A^g + Kg] u, » 0 (3) 

[AgiJ Agg - A^ I] U| * 0. (4) 


This work was supported by the NASA-Langley Research Center under grant number NSG-1519. 
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for a given eigenvalue asslgiwcnU U is thus necessary 
that each lie tn the null space of 

Nj • CAj,! Agj . X, 0 (5) 



for 1 ■ li 2i n as shown by equation (4), Equa- 
tion (3) Is then regarded as a constraint on K. once u. 
has been determined. If an independent sat of u.'s 
satisfying equation (4) is foundi thon U • ' 

CU|i Ug. •••i Ujj] will bo a modal matrl* satisfying 

equation (2), 

Lot R| be a matrix whose colu;ans form an ortho- 

normal basis for the null space of fl^. Then the 

columns of all possible modal matrices satisfying 
equation (2) can be expressed as linear ccmblnations 
of the columns of all Rj's. Let P denote the desired, 

although not necessarily realizable, modal matrix given 
as a design specification. The eigenvector selection 
Is made by projecting P, column by column, onto the 
vector space containing all possible modal matrices. 
This results in n sets of constants, a. * 

[oj^, 02 ^. < 1*1,2, •*•, n such that 

• Rj Op If the resulting eigenvectors U| are 

linearly independent, then the realizable modal matrix 
U which Is closest to P In a Euclidean norm sense has 
been obtained. If the resulting eigenvectors are not 
linearly Independent or If the resulting modal matrix 
Is unsatisfactory, the ojj parameters are perturbed 

slightly by weighting the projection norm. Once an 
Invertible modal matrix Is obtained, the closed-loop 

system matrix Is computed as A • UAU’\ The feedback 
gain matrix K Is then calculated from 

K - {[A„! A^g] • [An • Ai23) (6) 

where A Is partitioned to gonform to A. 


Suppose now that the elements of the plant 
matrices A and B are dependent on a scalar parameter 
p. Then the closed-loop system Is also dependent on p 
and the extent of this dependency, that Is the sensi- 
tivity of the system eigenvalues and eigenvectors with 
respect to p, Is related to the choice of modal matrix 
U. The second stage of the elgensystom assignment 
problem consists of reducing the elgensystem sensi- 
tivity by appropriate modification of the initial 

modal matrix. Let ^ and ^ denote the sensitivities 

of the plant parameters with respect to p. Then 


^ , dA . ^ . 
dp dp dp 


(7) 


and the eigenvalue sensitivities are given by [4] 



O “i 


while the eigenvector sensitivities are 


( 8 ) 


du^ 


z 

j*l 


f 


1 

_ ^i ' ^j _ 



(9) 


where V ■ [v^, Vg, •••, v^] Is a normalized reciprocal 
basis for A. A sensitivity cost function is defined by 


where tj, Wj, 1*1,2, ••• , n are positive weighting 

constants. Modifications to the initial modal matrix 
are calculated using a gradient search proceduru in 
which the sensitivity cost function d is reduced by 
choice of the constants Ojj. The sensitivity reduction 

stage of the design process is highly designer inter- 
active. At tho completion of each iteration a deci- 
sion to continue or terminate the sensitivity reduc- 
tion must be made. This decision often Involves a 
compromise between achieving low elgensystem sensi- 
tivity and retaining desired modal behavior. 


in. THE SERVO PROBLEM 


The procedure outlined above for the elgensystem 
assignment problem is applicable to output feedback 
and servo problems if observers are utilized to 
realize the feedback law obtained under the assumption 
of complete state accessibility. The design procedure 
for the servo problem Is developed here. Let the 
plant be described by it ■ Ax + 3u, y » Cx where 

X c r”, u e r"’, y e R*". In addition, the plant Is 
assumed to be completely state controllable and ob- 
servable, For clarity o^ presentation In this section 
It Is assumed that a change of coordinates has been 
made such that C Is In rank reduced form, C ■ [Ij,: 0], 

The block diagram In Flguro 1 Illustrates the 
Integral plus proportional feedback control structure 
to be considered, The dotted line In this figure in- 
dicates what the feedback structure would be If all 
plant states were available for measurement. The first 
step In the servo problem design procedure is to select 
gain matrices and Kg to achieve desirable system 

performance. The state dynamics of the system In Figure 
1 are given by 


A + BK, 
-C 


BK- 


Cu]. 


(n) 


It Is easily seen that the state matrix of this plant/in- 
tegrator system can be written as 


‘ A 1 0“ 
1 

X 

'b' 

.-c i o_ 

T 

.0_ 


CK. i Kg]. 


( 12 ) 


Since this matrix has the form of equation (1), the pro- 
cedures of the elgensystem assignment problem are appli- 
cable to the determination of K, and K, if the rank of 
the matrix ' ^ 



AB 

a^b| 

...j a""'"1b‘ 

_0 

CB 

cab! 

... 1 ca'’^'”^b_ 


is n+r. If this rank condition is not met, adjustment 
of the number of inputs and/or outputs may correct the 
situation. Therefore, given a designer specified 
eigenvalue matrix Ap and desired modal matrix Pp, the 

elgensystem assignment procedure yields gain matrices 
K^ and Kg which realize the specified eigenvalue matrix 
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Ap ind ) ipfldal iwtrfx Up» As noted fn Section II» Up 

usually reuresents a designer selected caiprotnlse 
between rial liable nsdal matrix nearest Pp and the 

real liable modal matrix resulting In Htinimuiti elgen- 
systoin sensitivity local to Pp. 

The second step in the servo prcbleiii design pro- 
cedure Is the construction of a reduced order observer 
to estimate those states which are not accessible at 
the output of the plant. A block diagram of a reduced 
order observer Is shown within dottod lines in Figure 
2. The observer system is described by £5] 

i ■ Fa ♦ Gy ♦ TOu (13) 

where 

^ * ^22 *^^12* ()^) 

T-C-Lil], (15) 

Q “ (A22 * ^ ^12^^ ^ ^^21 * ^ ^11^ 

,..’d A is »» '' plant matrix partitioned to conform with 
t.i wi C as 


[^21 '/ 22 J 

Furthermore, the transpose of equation (14), 

fT . <45 - aJ lT , (U) 


has the form of equation (10) so that the eigensystem 
assignment procedure may be utilized to determine L 
given desired eigenvalue and eigenvector matrices, 

Ao ’’o' 

The final step In the servo problem design pro- 
cedure Is the determination of the feedback gain 
matrices R and Kj in the complete integral plus pro- 
portional servo control system of Figure 2. This 
system Is described by the equation 


«• • 


- 

1 

1 


“ - 


- - 

• 

X 


A+BRC 

|.!^2. 

1 8K3 


X 


0 

• 

a 

-C 

! 0 

1 0 


’‘c 


I 

V 



T 

T ■ ^ - 


w 



• 

z 


GC+T8RC 

•|TBK2 

1 2+TBK3 


Z 


0 


(18) 


Then fore, choosing 


the state matrix of 


It follows that the 


K, ■ fiC 1 KgT , 
equation (20) becomes 


A « 


• A I f 

Ap I Af. 

Lo :f . 


eigenvalues of A are given by 

I 


A • 


A ; 0 

mmwttm 

,0 »Ap, 


( 21 ) 

(22) 

(23) 


and, assuming that the observer eigenvalues arc chosen 
distinctly from the eigenvalues of the plint/lntegrator 
system, the eigenvectors are given by 


U - 


U„ i QU, 


0 I U„ 


(24) 


In the original coordinates, the system modal matrix is 


U ■ 


U 

RU. 


B. 


Wc 

su. 


(25) 


where Q, R, S are linear transformation matrices. Thus, 
the servo system of Figure 2 realizes an eigensystem in 
which the plant/integrator dynamics and observer dynamics 
are decoupled and can be independently specified. Par- 
titioning K, of equation (12) as 


Kj ” [Kjjl > ( 26 ) 


where K^| consists of the first m colum.ns of K|, the 
solution for Kj and R from equation (19) is given by 

K 3 • K ,2 

and 

R ■ + K,2 L • (28) 

A computer program based on this design procedure 
has been written and is currently being evaluated. 
Initial results indicate that servo control systems 
having satisfactory modal behavior and low eigensystem 
sensitivity to changes in plant parameters can be ef- 
fectively designed by this procedure. A difficulty of 
this procedure, and a subject of future research, is 
the uncertainty of how to translate system performance 
specifications into specification of a desired modal 
matrix. A tutorial example Is presented In the next 
section to illustrate the new design procedure. 


Making the coordinate transformation 


’ X ’ 


1 

0 

0 


X 

*c 

■ 

0 I 0 


"c 

_ w _ 


.-T 0 I. 


.a . 


equation (18) becomes 


• 

X 


A+B(RC+K3T) 

BKj ; BK3 


X 


• 

■ 

-C 

0 1 0 



+ 

A 

■■ ■ 


0 

0 1 fJ 


V> 

w 

m ■■ 



(19) 
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IV. A DESIGN EXAMPLE 


Consider the system x » Ax + Bu, / » Cx with 
X c R'’, u c R^ and y c R^ where 


A ■ 


■- 2.0 

0.0 

1.0 


0.0 1.0 

-(2.0)/p (1.0)/p , 

1.0 -2,0 . 


B » 


I 

0 


( 20 ) 


1.0 

0.0 

0,0 


0.0 

(1.0)/p 

0.0 : 


C » 

1.0 

0.0 

0.0 


. 0.0 

1.0 

0 . 0 , 


and p ts a plant parameter with nominal value 1.0. It 
Is desired to design a servo control system so that 

T 

y “ [y] ^iV tracks step reference Input signals with 
zero steady-state error. In addition, the outputs y^ 
and yg should be decoupled, have Httle or no over- 
shoot, and have dominant time constants of approxi- 
mately 1 second and 1/2 second, respectively. 

In order to meet the desired tracking requirements 
for step reference Inputs, It Is decided to use a servo 
control system Incorporating I'xteqral error feedback as 
shown 1n Figure 2. A reduces order observer of order 
one Is Incorporated to estimate state which Is not 

available at the output of the plant. To meet the re- 
sponse time constant specifications, the closed loop 
eigenvalues are assigned as 

■ -1.0 Xg » -2.0, X3 » -8.0 

■ -9.0 Xg * -10.0, -7.0 

_ Finally, a diagonally dominant modal matrix specifica- 
’ tion Is made to provide decoupling between the out- 
puts. 

The design procedure consists of first calcula- 
ting the feedback matrices and K, of Figure 1 

assuming complete state feedback. A reduced order 
observer Is then constructed to estimate those states 
not actually available at the output. Finally, feed- 
back gain matrices K 3 and R are computed. The outcome 

of this design procedure Is given below. 


F “ 

[ 

-7.0] 


TB » 

c 

-5.0 

0.0] 

G - 

c- 

•24.0 

1.0] 

R » 

t 

•152.0 

33.2 

-5.4l 

-8.8j 


r 

39.1 

20. 9 I 

h 

L 

-4.8 

15. 9j 


The complete system Is described by the equation x • 
Ax + Bu where 



r - 154.0 

-5.4 

1.0 

39.1 

20.9 

- 27.4 


33.2 

- 10,8 

1.0 

- 4.8 

15.9 

6.3 

A 

A 

1.0 

1.0 

- 2,0 

0.0 

0.0 

0.0 

A “ 

‘ - 1.0 

0.0 

0.0 

0,0 

0.0 

0.0 


0.0 

- 1.0 

0.0 

0.0 

0.0 

0.0 


L 736.0 

28.1 

0.0 

- 195.5 

- 104,5 

129.9 


and 

r 0.0 0.0 *1 

0,0 0.0 
S , 0.0 0.0 

“ 1.0 0.0 ‘ 

0.0 1.0 

L 0.0 0,0 J 

The system response for Input reference signals 


"1 ' 

A 

’1 ■ 


and u » 


.0. 


.1 - 


are shown In Figures 3a and 3b, respectively. 


V. REFERENCE.S 


[1] R. A. Maynard, R. R. Mielke, S. R. Liberty, 
and S. Srlnathkumar, "A New Spectral Synthesis Pro- 
cedure for Multivariable Regulators", Proc., Allerton 
Conference on Communication, Control, and Computing, 
pp, 754-761, October 4-6, 1978. 

[2] R. R. Mielke, R. A. Maynard, S. R. Liberty, 
and S. Srinathkuirar, "Design of Robust State Feedback 
Controllers Via Eigenvalue/Eigenvector Assignment", 
Conference Record, Twelfth Asilomar Conference on Clr^ 
cults, Systems, and Computers, pp. 287-291, November 
6 - 8 , 1979. 

[3] S. R. Liberty, R. R. MieHe, and R. A. May- 
nard, Final Report on NASA-Langley Research Center 
grant number NSG-1519, in preparation. 

[4] J. H. Wilkinson, The Algebraic Einenvalue 
Problem , London: Oxford University Press, 1965. 

[5] 0. J. Luenberger, "An Introduction To Ob- 
servers", IEEE Transactions on Automatic Control, 

V 0 I. AC-16, No. 6 , pp. 596-602, December 1971. 



• FIGURE 1 

SERVO CONTROL SYSTEM WITH COMPLETE STATE FEEDBACK 
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